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Abstract 

Let So be a smooth and compact real variety given by a reduced 
regular sequence of polynomials fx,...,f p . This paper is devoted to 
the algorithmic problem of finding efficiently a representative point for 
each connected component of Sq . For this purpose we exhibit explicit 
polynomial equations that describe the generic polar varieties of So ■ 
This leads to a procedure which solves our algorithmic problem in time 
that is polynomial in the (extrinsic) description length of the input 
equations fx, • • ■ , fp and in a suitably introduced, intrinsic geometric 
parameter, called the degree of the real interpretation of the given 
equation system fx, . ■ ., f p . 
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1 Introduction 

The core of this paper consists in the exhibition of a system of canoni- 
cal equations which describe locally the generic polar varieties of a given 
semialgebraic complete intersection manifold So contained in the real n- 
dimensional affine space 1R™ . This purely mathematical description of the 
polar varieties allows the design of a new type of efficient algorithm (with 
intrinsic complexity bounds), which computes, in case that So is smooth 
and compact, at least one representative point for each connected compo- 
nent of So (the algorithm returns each such point in a suitable symbolic 
codification). This new algorithm (and, in particular, its complexity) is the 
main practical outcome of the present paper. Let us now briefly describe 
our results. 

Suppose that the real variety So is compact and given by polynomial equa- 
tions of the following form: 

fi(Xi, . . . , X n ) = ■ ■ ■ = f p (Xi, . . . , X n ) = 0, 

where p, n G IN, p < n and fi,...,f p belong to the polynomial ring 
Q [Xi , . . . , X n ] in the indeterminates X\ , . . . , X n over the rational num- 
bers Q . Let d be a given natural number and assume that for 1 < k < p 
the total degree deg/^ of the polynomial is bounded by d. More- 
over, we suppose that the polynomials fi, ■ ■ ■ , f p form a regular sequence in 
Q [Xi , . . . , X n ] and that they are given by a division- free arithmetic circuit 
of size L that evaluates them in any given point of the real (or complex) 
n -dimensional affine space R n (or C n ). Further, we assume that the Ja- 
cobian J(/i, . . . , f p ) of the equation system f\ = ■ ■ ■ = f p = has maximal 
rank in any point of So (thus, implicitly, we assume that So is smooth). 
Let Wo := V(f\, . . . , f p ) denote the (complex) algebraic variety defined by 
the polynomials fi,...,f p in the affine space C™ . We denote the singular 
locus of Wq by SingWo ■ 

Moreover, let us suppose that the variables X±, . . . ,X n are in generic po- 
sition with respect to the equation system fi,...,f p . For 1 < i < n — p 
let Wi be the i -th formal (complex) polar variety associated with Wo (and 
the variables X p+ i, . . . , X n ). 

Further, let us denote the real counterpart of Wi by Sj := WinlR™ . We call 
Si the i— th formal real polar variety associated with the real semialgebraic 
variety So (and the variables X p+ i, . . . , X n ) . It turns out that the (locally) 
closed sets Wi\SingWo (resp. Sj ) are either empty or complex (resp. real) 
smooth manifolds of dimension n — (p + i) . Moreover, for 1 < i < n — p , 
one sees easily that 

Wi := Wi \ SingWo 

is the z— th polar variety (in the usual sense) associated with Wo and the 
variables X p+ i,...,X n (here, Wi \ SingWo denotes the Q-Zariski closure 
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in (D n of the quasi-affine variety Wi \ SingWo ). For a precise definition of 
the notion of formal polar varieties and of polar varieties in the usual sense 
we refer to Section 2. 

Suppose that the real variety Sq is non-empty and satisfies our assumptions. 



In Theorem 10 of this paper we show that every real polar variety S{ = 
Wi flE", 1 < i < n — p, is a non-empty, smooth manifold of dimension 
n — p — i containing at least one point of each connected component of the 
real variety Sq . In particular, the real variety S n - P is a finite set containing 
at least one representative point of each connected component of So . 
Under the same assumptions we show in Theorem || that for 1 < i < 
n—p the quasi-afnne variety Wj\ SingWo is a locally complete intersection 
that satisfies the Jacobian criterion. More precisely, the quasi-afnne variety 
Wi\SingWo is a smooth manifold of codimension p+i that can be described 
locally by certain regular sequences consisting of the polynomials fi,---,f p 
and i many well-determined p-minors of the Jacobian J(fi, ■ ■ ■ , f p ) of 
the fi,...,f p . In particular, the quasi-afnne variety W n -p \ SingWo is 
zero-dimensional, whence W„_ p = W n - P \ SingWo ■ Thus W n - P is a zero- 
dimensional complex variety that contains a representative point of each 
connected component of the real variety So . 

The practical outcome of Theorem || and Theorem |l(] consists in the design 
of an efficient algorithm (with intrinsic complexity bounds), which adapts 
the elimination procedure for complex algebraic varieties developed in (3(J 
and |3lJ to the real case. Under the additional assumption that for any 
1 < k < p, the intermediate ideal (/i, •••,/&) generated by fi,...,fk 
in Q [Xi,...,X n ] is radical, we shall apply this procedure to the p(^L^) 
well-determined equation systems of Theorem ||, which describe the zero- 
dimensional algebraic variety W n - P = W n _ p \ SingWo locally. In order to 
find at least one representative point for every connected component of the 
real variety Sq , we have just to run the procedure of []30| and on all 
these equation systems. Counting arithmetic operations in Q at unit costs, 
this can be done in sequential time 



\L(ndS) 



O(i) 



/ n 

where 5 is the following geometric invariant of the regular sequence /i , 



ip ■ 



5 := max{max{degU(/i, ... ,f k )\ SingW \l <k<p}, 
maxjdeg Wi|l < i < n — p}} 



(here, deg V(f±, . . . , fj~) \ Wo and degWj denote the geometric degree in 
the sense of p7| of the corresponding algebraic varieties). 



This is the content of Theorem |ll| below. For any 1 < k < p and any 



1 < i < n — p the quantity 5 bounds the degree of the algebraic variety 
V(fi, . . . , fk) \ SingWo and of the i-th. polar variety Wi = Wi\ SingWo ■ 
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In |30| and j3l| the quantity maxjdeg V(f±, . . . , fi)\l < i < p} is called 
the geometric degree (of the complex interpretation) of the equation system 
fi,...,f p . In analogy to this terminology, we shall call 5 the geometric de- 
gree of the real interpretation of the equation system fx , . . . , f p . In view of 
the complexity result above we shall understand the parameter 5 as an in- 
trinsic measure for the size of the real interpretation of the given polynomial 
equation system. Nevertheless, the word "intrinsic" should be interpreted 
with some caution in this context: observe that the complexity parameter 
5 depends rather on the equations fi, ■ ■ ■ , f p and their order than just on 
the variety Wo \ SingWo . 

In order to make our complexity result more transparent we are going now to 
exhibit, in terms of extrinsic parameters, some estimations for the intrinsic 
system degree 5 . 

Let us write d\ := deg/i, . . . , d p := deg/ p and let D := d\- ■ ■ d p denote 
the classical Bezout number of the polynomial system fi,...,f p . Then 
we have the following degree estimations for the complex algebraic variety 
W = V(f 1 ,...,f p ) 

deg S < deg W < D < a? 

(So denotes again the Q-Zariski closure in C n of the real variety So )• 
On the other hand, we conclude from Theorem || that, for every i, 1 < 
i < n — p, the polar variety Wi is defined by the initial system fi,...,f p 
and certain p -minors of the Jacobian J(/i, . . . , f p ) ■ Let us denote the 
maximum degree of these p -minors by Cj . It turns out that for, any 1 < % < 
n — p, the polar variety Wi is a codimension one subvariety of Wj_i . Now 
one sees easily that the quantity D{ := D c\ ■ ■ ■ Cj represents a reasonable 
"Bezout number" of the variety Wj and that this Bezout number satisfies 
the estimate deg W« < Di . Putting all this together, we deduce the following 
estimate for the intrinsic system degree 5 : 

5 < T>n- V = Dei ■ ■■c n - p . 

Observing that for any i, 1 < i < n—p, the inequality Cj < d\ + • • • + d p — p 
holds, we find the estimations: 

5 < D(di + --- + d p - p) n ~ p < a?{pd - p) n ~ p < p n - p d n . 

In conclusion, our new real algorithm has a time complexity that is, in 
worst case, polynomial in the "Bezout number" Dc\ ■ ■ ■ c n - p of the zero- 
dimensional polar variety W„_ p . 

Our complexity bound ( "]) L(nd5)°^' depends on the intrinsic (geometric, 
semantic) parameter 5 and on the extrinsic (algebraic) parameters d and n 
in a polynomial manner, and it depends on the syntactic parameter L only 
linearly. In this sense one may consider our complexity bound as intrinsic. 
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Our real algorithm promises therefore to be practically applicable to special 
equation systems with low value for the intrinsic parameter 5 . 

On the other hand, even in worst case our algorithm improves upon the 
known dP^ -time procedures for the algorithmic problem under consider- 
ation, also in their most efficient versions [||, || (see also ||, [p^| , fl9| , 3£], 
Q, Q, @, 0). However, this distinction does not become 

apparent when we measure complexities simply in terms of d and n (all 
mentioned algorithms have worst-case complexities of type d°^ ), but it 
becomes clearly visible when we use the "Bezout number" just introduced as 
complexity parameter. Only our new algorithm is polynomial in this quan- 
tity. On the other hand, we are only able to reach our goal of algorithmic 
efficiency by means of a strict limitation to a purely geometric point of view. 
For the moment there is no hope that any of the standard questions of real 
algebra (e.g. finding generators for the real radical of a polynomial ideal or 
the formulation of an effective real Nullstellensatz) can be solved within the 
complexity framework of this paper (compare |52| , j7j and |8| ) . 

In conclusion we may say that this paper establishes a new connection be- 
tween the algorithmic complexity of finding a representative set of real solu- 
tions of a given polynomial equation system and the geometry of the (com- 
plex) algebraic variety defined by this system. However, there is a price to 
pay for that: this connection becomes only visible if we restrict ourselves 
to reduced complete intersection systems that define smooth, compact real 
varieties. 

Our (algorithmic and mathematical) methods and results represent a non- 
obvious generalization of the main outcome of jjj], where an intrinsic type 
algorithm was designed for the problem of finding at least one representative 
point in each connected component of a real, compact hypersurface given by 
an n-variate, smooth polynomial equation / of degree d>2 with rational 
coefficients (such that / represents a regular equation of that hypersurface). 
This is the particular case of codimension p = 1 of the present paper, and 
our setting leads to the complexity bound L^dS) ^ proved in ||]. 



2 Polar Varieties 

2.1 Notations, Notions and General Assumptions 

Let X\, . . . ,X n be indeterminates (or variables) over the rational numbers 
Q and let polynomials fi,---,f p € Q [Xi, . . . , X n ] with 1 < p < n be 
given. Let (D n and IR n denote the n -dimensional affine space over the 
complex and the real numbers, respectively. We think C n to be equipped 
with the Q -Zariski topology, whereas, on JR n , we consider the strong 
topology. For any subset U C C n we denote its Q -Zariski-closure by 
U . By X := (Xi, . . . , X n ) we denote the vector of variables Xi, . . . , X n 
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and by x := (x\, . . . ,x n ) any point of the affine space C n or R™. We 
suppose that the polynomials fi, ■ ■ ■ ,f p form a reduced regular sequence in 
Q [Xi , . . . , X n ] (here " reduced" means that for any 1 < k < p the ideal 
(/i, . . . , fk) is radical). The Jacobian of these polynomials is denoted by 



For any point x € € n we write 



dXj 



dfk 
dXi 



l<k<p 
l<j<n 



(x) 



l<k<p 
l<j<n 



for the Jacobian of the polynomials /i , . . . , f p at x . 

The common complex zeros of the polynomials f±, . . . , f p form an affine, 
Q-definable subvariety of C n , which we denote by 



W n 



V(h,...,f p ) :={xeV n \h{x) 



x 



0}. 



A point x G Wo = V(fi, . . . , f p ) is said to be non-singular or smooth (in 
Wo ) if the rank of the Jacobian of /i , . . . , f p in x is p . Otherwise x is 
called a singular point of Wo ■ By SingWo we denote the set of all singular 
points of Wo ■ Since we suppose that the ideal (f±, . . . , f p ) is radical, our 
notion of a smooth point coincides with the usual one for algebraic varieties 
by Jacobi's criterion. 

Remark 1 

If x G Wo is smooth, then the hypersurfaces defined by the polynomials 
fi , . . . , f p intersect transversally at the point x . 



Definition 2 

For every i, 1 < i < n — p, let Aj denote the set of all common complex 
zeros of all p -minors of the Jacobian J(/i, . . . , f p ) corresponding to the 
columns {1, . . . ,p + i — 1} . In other words, A» is the determinantal variety 
defined by all p-minors of the submatrix Jf +l_1 (/i, • • • , f p ) determined by 
the columns {1, . . . ,p + i — 1} of the Jacobian J(/i, . . . , f p ) ■ 
We introduce the affine variety 

Wi := Wo n A, 

associated with the linear subspace of (D n , namely 

X^- 1 := {x G € n \X p+i (x) = ... = X n (x) = 0} 

and call Wi the i -th formal polar variety of Wq ■ 

By _ 

Wi-.= Wi\ SingWo 

we denote the i-th polar variety (in the usual sense) of the variety Wo . 



2 POLAR VARIETIES 



7 



Remark 3 

• Our definition of polar and formal polar variety depends rather on the 
regular sequence fi,---,f p than on the algebraic variety Wq . The 
ad hoc term "formal polar variety" is only used in this paper for the 
purpose of clarification of our subsequent mathematical arguments. 

• The index i reflects the expected codimension of the polar variety Wi 
in Wo . With respect to the ambient space V 1 , the expected codimen- 
sion of Wi is p + i (see Theorem || below for a precise statement). 

• According to our notation, the common zeros of all p -minors of the 
Jacobian J(fi, ■ ■ ■ , f p ) form the determinantal variety A n _ p+ i . Ob- 
viously, we have Sing Wo = Wo PI A n _ p+ i = W n - p+ \ . 

• The formal polar varieties Wi, 1 < i < n — p, constitue a decreasing 
sequence. In particular, we have 

Wo D Wi D ■ ■ ■ D Wi D • • • D W n -p D W n - P+1 = Sing Wq. 

The concept of polar variety goes back to J.-V. Poncelet. Its development 
has a long history: Let us mention among others the contributions of F. 
Severi, J. A. Todd, S. Kleiman, R. Piene, D. T. Le, B. Teissier, J. -P. Henry, 



M. Merle ... (see e.g. [58] and the references quoted there). 



2.2 Local Description of the Determinantal Varieties 

In this subsection we develop a succinct local description of the determinan- 
tal varieties Aj, 1 < i < n — p. The following general Exchange Lemma 
will be our main tool for this description (this lemma is used in a similar 
form in |32j). It describes an exchange relation between certain minors of a 
given matrix . 

Let A be a given (p x n) -matrix with entries a^- from an arbitrary com- 
mutative ring. Let I and k be any natural numbers with / < n and 
k < min{p, 1} . Furthermore, let Ik := (h, ■ ■ ■ , it) be an ordered sequence of 
k different elements from the finite set of natural numbers {1, . . . , 1} and let 
Ma (Ik) '■= Ma(£1) ■ ■ ■ ,ik) denote the A; -minor of the matrix A built up by 
the first k rows and the columns i\, . . . , if- . If it is clear by the context what 
is the matrix A , we shall just write M(i\, . . . , i^j := M {Ik) '■= Ma (Ik) ■ 

Lemma 4 (Exchange Lemma) 

As before let a matrix A and natural numbers I and k be given, as well 
as two intersecting index sets Ik = (h, ■ ■ ■ , ik) and Ik-i = (ji, ■ ■ ■ ,jk-i) ■ 
Then, for suitable numbers £j E {1,-1} with j € Ik\ Ik— I we have the 
following identity: 

(*) M(4_ 1 )M(J fe )= £ £j M(I k \{j})M(I k ^U{j}). 
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Proof 

Consider the following ((2k — 1) x (2k — 1)) -matrix L with entries from the 
given matrix A: 





Li(h) 









Lk-i(h) 


Li(h-i) 


Li(h) 


Lk(Ik-i) 


Lk(h) 



Here, for any 1 < j < k, Lj (I^) denotes the row vector of length k that 
we obtain selecting, from the j -th row of the matrix A , the k elements 
placed in the columns Ik = (h, ■ ■ ■ ,ik) ■ Similarly, Lj (Ik-i) is obtained 
from the j -th row of A selecting the k — 1 elements placed in the columns 

4-i = (ii, • • • ,jk-i) ■ 

Now it is not difficult to verify the identity (*) by calculating the determi- 
nant det L of the quadratic matrix L via Laplace expansion in two different 
ways. First, by expansion of det L according to the first k — 1 columns of 
L , we obtain the left-hand side of (*) , disregarding the sign. Expansion of 
det L according to the first k — 1 rows of L leads to the right-hand side of 
(*) . This implies the identity (*) for an appropriate choice of the signs Ej , 
with j G 4 \ 4-1 • D 

Let m € Q [X±, . . . , X n ] denote the (p — 1) -minor of the Jacobian 
J(/i, . . . ,f p ) given by the first (p — 1) rows and columns, i.e., let 



m := det 



dfk 
dXj 



i<*:<p-i 
i<j< P -i 



We consider the determinantal variety A« outside of the hypersurface 

V(m) := {xe€ n \ m(x) = 0} 
and denote this localization by (Aj) m , i.e., we set 

(Ai) m :=Ai\V(m). 
From now on , for 1 < i\ < ■ ■ ■ < i p < n , let us denote by 

M(h, ...,i p ) 
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the polynomial in Q [X±, . . . ,X n ] denned as the p -minor of the Jacobian 
J(/i, . . . , f p ) built up by its p rows and the columns i±, . . . , i p . As before, 
we denote by 

M(h, . . .,i p )(x) 
the specialization of M(i±, . . . ,i p ) in a given point x £ <C n . 

Proposition 5 

Let 1 < i < n — p be arbitrarily fixed, and let m be the (p — 1)- minor 
defined above. Then the determinantal variety A, is locally (i.e., outside of 
the hypersurface V(m) ) described by the i polynomials 

M(l, ...,p- l,p), M(l, . . . ,p - l,p + 1), . . . , M(l, . . . ,p - l,p + i - 1). 

In other words, we have 

(Ai) m := {x G (TJ 1 ! m(x) + 0, M(l, . . . ,p-l, s)(x) = 0, s £ {p, . . . ,p+i-l}}, 

where M(l, . . . ,p — 1, s) denotes, as above, the p-minor of the Jacobian 
J(/i, . . . , / p ) buiJt up by the first p — 1 columns and the s -th column. 

Proof 

It suffices to show that 

(A t ) m D{i£ <D n | m(a;) ^ 0, M(l, . . . ,p — l,s) = 0,s €. {p, . . . ,p + i — 1}} 
holds. 

Let x* S (D n be any point satisfying the conditions m(x*) ^ and 
M(l, ... ,p — 1, = for every s G {p, . . . ,p + i — 1} . We have to 

verify that 

M(i 1 ,...,i p )(x*) = 

holds for all ordered p -tuples . . . , i p ) of elements of {1, . . . ,p + i — 1} . 
Applying the Exchange Lemma to m = M(l, . . . ,p— 1) and M(i\, . . . ,i p ) , 
we deduce the identity 

m(x*)M{i 1 ,...,i p ){x*) = 
e, M({i 1 ,...,i p }\{j})(x*)M(l,...,p-l, 3 )(x*) 

j'e{ii,...,i p }\{l,...,p-l} 

for suitable numbers ej G {1, —1} with j G . . . , i p } \ {1, . . . ,p — 1} . By 
assumption we have m(x*) 7^ and M(l, . . . ,p — l,j)(x*) = for all j G 
{p, . . . ,p + i — 1} . This implies that x* belongs to the set (Aj) m . 

□ 
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Notation 6 

In the sequel we shall simply write Mj for the p- minor M(l, . . . ,p — 
l,j) given by the first p—1 columns of J(/i,...,/ p ) and the column 
j e{p,...,n} . 

Remark 7 

• Proposition implies that the codimension of Aj outside of the hy- 
persurface V(m) is at most i . 

• Proposition [| holds also for the determinantal variety A n _ p+ i that 
defines the singular locus Sing Wo = W n - p+ i of the variety Wq . 
Hence, for any point x* £ (D n satisfying the condition m{x*) ^ and 
the n — p + 1 equations 



the Jacobian J(/i, . . . , f p )(x*) becomes singular. 

• Replacing the previously chosen (p—1) -minor m by any other (p—1)- 
minor of the Jacobian J(/i, . . . , f p ) , the statement of Proposition |^ 
remains true mutatis mutandis. 

2.3 Local Description of the Formal Polar Varieties 

The aim of this subsection is to show the following fact: 

Let the variables Xi,...,X n be in generic position with respect to the 
polynomials f\, . . . , f p , and let m be any (p — 1) -minor of the Jacobian 
J(/i, . . . , f p ) . In this subsection we are going to show that any formal polar 
variety Wi, 1 < i < n — p, is a smooth complete intersection variety outside 
of the closed set SingWo U V(rh) . Moreover, we shall exhibit a reduced 
regular sequence describing this variety outside of SingWo U V(m) . 

As in the previous subsection, let m £ Q [X%, . . . , X n ] denote the (p — 1) - 
minor of the Jacobian J(/i, . . . , f p ) built up by the first (p — 1) rows and 
columns. 

Let Y\, . . . , Y n be new variables and let Y := (Yi, . . . , Y n ) . For any linear 
coordinate transformation X = AY , with A being a regular (ra X n) — 
matrix, we define the polynomials 



M j (x*) = 0, j € {p,...,n}, 



Gi(r) := h(AY), G P {Y) := f p (AY). 



The Jacobian of G%, . . . , G p has the form 



J(G%, . . . , G p ) : 



dG k 



J(fi,...J P )A. 




l<k<p 
l<j<n 
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Using a similar notation as before, we denote by 

M(ix, ...,i p ) 

the p-minor of the new Jacobian J(G\, . . . ,G P ) that corresponds to the 
columns 1 < i\ < ■ ■ ■ < i p < n . 

Moreover, we denote by Mj the p- minor M(l, . . . ,p — determined 
by the fixed first p — 1 columns of J(G±, . . . ,G P ) and the column j G 
{p,...,n}. 

For p < r,t < n let Z r j be a new indeterminate. Using the following 
regular (n — p + 1) x (n— p+1) -parameter matrix 



Z :-- 



1 





1 















3 



p+i-l,p -^p+j-ljp+l 



->p+i,p+l 



J n,p+1 



->p+i,p+i-l 



J n,p+i—l 



J n,p+i ^n,p+i+l 



we construct an (n x n) -coordinate transformation matrix A := A(Z) , 
which will enable us to prove the statement at the beginning of this subsec- 
tion. 

For the moment, let us fix an index 1 < i < n — p . We consider the formal 
polar variety Wj outside of the hypersurface V(m) . Corresponding to our 
choice of i , the matrix Z may be subdivided into submatrices as follows: 



Z 



Oi,n—p- 



i+l 



Here the matrix is defined as 



Z 



p+i,p 



J np 



->p+i,p+i-l 



J n,p+i—l 



(i) (i) 

and Z\ and Z\ denote the quadratic lower triangular matrices bordering 
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2fW in Z , and Oj )n _ p _j + i is the i x (n — p — i + 1) zero matrix. Let 



A := A(Z) := 



Ip-i 
Oi :P -i 

On—p—i+l,p—l 



Op—\,i O p—\,n—p—i-\-\ 



Oi.n—v— 



n—p~i+l 



Here the submatrices I r and O r)S are unit or zero matrices, respectively, 



of corresponding size, and Z®,Z[ %) , and Z%' are the submatrices of the 
parameter matrix Z introduced before. Thus, A is a regular, parameter 
dependent (n x n) -coordinate transformation matrix. 

Like the matrix Z , the matrix A contains 

(n — p)(n — p+l) 



parameters Z r ^ which we may specialize into any point z of the affine 
space <D S . For such a point z € C s we denote the corresponding specialized 



matrices by A(z) , Z{'(z), Z$'(z) and ZW(z) . 

We consider now the coordinate transformation given by X = AY with 
A = A{Z) and calculate the Jacobian J(G±, . . . ,G p ) with respect to the 
new polynomials Gi, . . . ,G P . Recall that the coordinate transformation 
matrix A depends on our previous choice of the index 1 < i < n — p . 

According to the structure of the coordinate transformation matrix A = 
A(Z) we subdivide the Jacobian J(/i, . . . , f p ) into three submatrices 



J(fi,---,fp) 



U V w 



with 



U :-- 



df k 



, V :- 



l<k<p 
l<j<p-l 



, W:~- 



l<k<p 



dfk 
dXi 



l<k<p 
p+i<j<n 



From the identity J(G\, . . . , G p ) = J(/i, ...,f p ) A we deduce that our new 
Jacobian is of the form: 



J(Gi, ... , G p ) 



dG k 



i<k< P 

l<j<n 



u vzf> + wz^ wz. 



(0 



We are interested in a local description of the i-th formal polar variety 
Wi = WqDAi outside of the hypersurface V(m) , where m is the fixed upper 
left (p — 1)- minor of the Jacobian J(/i, . . . , f p ) (and also of its submatrix 
U ) . Since the coordinate transformation X = AY leaves the submatrix U 
unchanged, the (p — 1) -minor m remains fixed under this transformation. 
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From Proposition ^ we know that the localized determinantal variety (Aj 
is described by the i equations 



0, . . . , M p+l _i = 0, 



and by the condition m ^ 0. The p- minors M p , . . . , M p +i-% defin- 
ing these equations are built up by the submatrix [ U V] of the Jaco- 
bian J(/i, . . . , f p ) . Under the coordinate transformation A{Z) the matrix 
[U V] is changed into the submatrix 



of the Jacobian J(Gi, 
into the p— minors 



u vzf ] + WZ W 

, G p ) and the p-minors M p , . . . , M p+ j_i 



are changed 



M, 



of the matrix 



u vzf + WZ® 



(**) 



M„,...,M, 



p+i—l 



This implies the matrix identity 



[M p , Mp+i-i] Zf + [M p+i , ...,M n ] Z®. 



For the previously chosen index 1 < % < n—p , the coordinate transformation 
X = A{Z)Y induces the following morphism of affine spaces: 



$i : C"xC s ^C p x (D\ 



defined by 

(x,z) H- 



$i(x, z) := I . . . , /p(x), M p (ac, z), . . . , M p+i _i(>, z) 



Consider an arbitrary point z € <D S . We denote by A? the determinantal 
subvariety of C n defined by all p -minors of the matrix 
U Vzf\z) + WZ {i \z)] (which is a submatrix of the new Jacobian ob- 
tained by specializing the coefficients of the polynomials G\, . . . ,G p into 
the point z € € s ). Writing Wf := W) PI Af , one sees immediately that the 
zero fiber $~ 1 (0) of the morphism <I>j contains the set 

(W[) m := W n (Af) m . 

In other words, for any arbitrarily chosen point z € C? , the zero fiber ^(O) 
of the morphism <£j contains the transformed formal polar variety Wf , 
localized in the hypersurface V(m) and expressed in the old coordinates. 
We are going now to analyze the rank of the Jacobian of the morphism $j in 
an arbitrary point (x, z) G (D™ x C with x € (Wf ) m . Using the subdivision 
of the parameter matrix Z into the parts Z-jp and zip , the Jacobian 

J(&i) of the morphism $j can be written symbolically as 



dX <9Z« 



dzf dZ, 



2 J 
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We have 



d$j d<z> t 

dX dZ® 



J(fli ■ ■ ■ i fp) Op yn —p—i+i 

dM p dM p 



dZ p+l , p !■■■■> dZ n 



0\ t n—p—i-\-\ 



Oi 



,n—p—i+l 



dM, 



p+i — 1 



dM, 



P+i— 1 



where the columns correspond to the partial derivatives of &i with respect 
to the variables 

Xi, . . . , X n , Zp+i^p, . . . , Z n ^p, . . . , Zp+i t p+i—±, . . . , Zn^p+i—i 

(in this order). The entries O r ^ denote here zero matrices of corresponding 
size and the row matrices labeled by " * " represent the partial derivatives 
with respect to the variables X±, . . . , X n of the minors M p , . . . , M p+ i-\ . 
These row matrices will be irrelevant for our considerations. 

d<S>, 



Furthermore, the third submatrix 



dM v 



dM p 

8Z p -\-\ tV ' ' ' ' ' dZ p +i-i, P 

Oi,*-i 

Oi,i-i 
Oi,i-i 



and the last submatrix 



dzf 

Op^i-2 
U -2 



8M P+1 



dZ, 



P+2,p+l 



of J($i) can be written as 




■ ■ ■ 



' dZ, 



p+i-l,p+l 



lti -2 
U -2 



9M p+i _i 





dZ, 



(0 



2 J 



of J(&i) is a zero matrix since the p- 



minors M p , . . . , M p+ j_i are indepedent of the parameters Z r t occurring in 



the submatrix Z 2 of the coordinate transformation matrix A(Z) . 

p (*i) is ( 



Therefore, the Jacobian J(3>i) is of full rank p + i wherever the submatrix 



dX dZ® 
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is of full rank p + i . On the other hand, considering the i row matrices 
contained in J($i) for p<j<p+i— 1 

dMj dMj ' 

_9Z p +i,j ' ' 9Z n jj ' 

we see that the representation (**) of the transformed p -minors Mj implies 
the identity 

dMj dMj 
dZp+ij ' ' dZ n j 

Thus, we obtain the representation 



Op,n—p— 
0\,n—p—i+\ 

I p+i ,...,M n ] _ 

Since all entries of the submatrix J(&i) of the Jacobian J(&i) belong to 
the polynomial ring Q [X\,. . . ,X n ], we see that the rank of the matrix 
J(<£j) in a given point (x,z) £ C" x C s with x G (Wi) m depends only on 
the choice of x . According to our localization outside of the hypersurface 
V(m) , let us consider an arbitrary smooth point x of Wq = V(/i, • • • , fp) 
satisfying the condition m(x) ^ 0. Suppose that the submatrix J(<E>j)(x) 
is not of full rank, i.e., that 

rk J($i)(x) < p + i 

holds. This latter inequality is valid if and only if all p -minors M p+i , . . . , M n 
of the Jacobian J(/i, . . . , f p ) vanish at x . Let z £ C s be any parameter 
point such that the pair (x,z) belongs to the fiber $~ 1 (0) of the mor- 
phism <£j . Since the p -minors M p , . . . , M p+ j_i of the transformed Jacobian 
J(Gi,...,G p ) must vanish at (x,z), we deduce from (**) that 

[0, . . . , 0] = [M p (x), Mp+i^x)} Z®{z) 

holds (here z[ l \z) denotes again the matrix obtained by specializing the 

(i) 

entries of Z\ into the corresponding coordinates of the point z G (D s ). 

(i) 

Because of the lower triangular form of the regular matrix Z\ , the latter 
matrix equation holds if and only if the conditions 

M v+i _ x {x) = ... = M p (x) = 0. 



= [M p+i ,...,M n }. 



<J(fii • • • i fp) Op yn —p—i+i 
* [M p+i ,...,M n 



O 



l,n— p— i+1 



[A 
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are satisfyied. Therefore, our assumptions on x and z imply m{x) ^ 
and Mp{x) = ■ ■ ■ = M n {x) = 0. However, by Remark [?] this means that the 
Jacobian J {fx, ■ ■ ■ , f p ){x) is singular. Hence, x is not a smooth point of 
Wo , i.e., x G Sing Wq , which contradicts our assumption on x . 

Now, suppose that we are given a point {x, z) G C n x C s that belongs to 
the fiber $~ 1 (0). Then x belongs to Wo- Further, suppose that x is a 
smooth point of Wo outside of the hypersurface V{m) . Let us consider the 
Zariski-open neighbourhood U of x, consisting of all points x G C n with 
m{x) 7^ and rk J{fx, • • • , f p ) = P, i-e., we consider 

U := (D n \ {Sing W U F(m)) . 

We are going to show that the restricted morphism 

is transversal to the origin 6 C x C . 

In order to see this, consider an arbitrary point {x, z) of U x CP that satisfies 
the equation $i{x,z) = 0. Thus, x belongs to U n Wq and is, therefore, 
a smooth point of Wo , which is outside of the hypersurface V{m) . By the 
preceding considerations on the rank of the Jacobian J{&i) it is clear that 
J(&i) has the maximal rank p + i at {x,z). This means that (x,z) is a 
regular point of <5j . Since (x, z) was an arbitrary point of $~ 1 (0)n (C/xC) , 
the claimed transversality has been shown. 

Now, applying the Weak-Transversality-Theorem of Thorn-Sard (see e.g. 
|P2| ) to the diagram 

*7 X (0) n (C7 x C s ) C n xC s 

\ I 

one concludes that there is a residual dense set fij of parameters z G s for 
which transversality holds. This implies that, for every fixed z 6 fl;, the 
transformed and localized formal polar variety 

Wi \ {Sing W UV{m)) 

is either empty or a smooth variety of codimension p + i . This variety can 
be described locally by the polynomials 

(* * *) h{X), f p {X), M P {X, 2O, ... , M p+i -x{X, z) 



that form a regular sequence outside of SingWo U V{m) . Up to now, our 
considerations concerned only the change of coordinates for an arbitrarily 
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fixed 1 < i < n — p. However, Q := Cii=i is a dense residual parameter 
set in C s from which we can choose a simultaneous change of coordinates 
for all 1 < i < n — p . For every choice z € and 1 < i < n — p the trans- 
formed formal polar variety Wf is, outside of the closed set SingWoUV (m) , 
a smooth complete intersection variety described by the (local) regular se- 
quence (* * *) . One sees now easily that the affine space H s contains a 
non-empty residual dense set of parameters z such that the conclusions 
above apply to the coordinate transformation X = A(z)Y . Moreover, z 
can be chosen from Q s . 

Taking into account Proposition || and Remark [7], we deduce the following 
result from our argumentation: 

Theorem 8 

Let Wo = V(fi, . . . , f p ) be a reduced complete intersection variety given 
by polynomials f\ , . . . , f p in Q [X\ , . . . , X n ] and suppose that the variables 
Xi, . . . , X n are in generic position with respect to /i, . . . , f p . Further, let m 
be the upper left (p — 1) -minor of the Jacobian J(/i, . . . , f p ) . Then, every 
formal polar variety Wi, 1 < i < n — p, localized with respect to the closed 
set SingWo U V(m) , is either empty or a smooth variety of codimension 
p + i that can be described by the equations 

hi-- - 1 fp,Mp, ■ ■ ■ ,M p+i -i, 

where Mj, p<j<p + i—l, is the p - minor of the Jacobian J {fx, . . . , f p ) 
given by the columns 1, . . . ,p — l,j . Then the polynomials 

/i) • • • ) fp, M p , . . . , M p+i _i 
form a regular sequence outside of SingWo U V(m) . 

Remark 9 

Taking into account that the argumentation on the localization with respect 
to the fixed (p — 1) -minor m remains valid mutatis mutandis for any other 
(p — 1) -minor rh of the Jacobian J {fx, . . . , f p ) , Theorem [| can be restated 
for any hxed {p — 1) -minor just by reordering of columns and rows of the 
Jacobian J (fx, • • • , fp) ■ 

2.4 Existence of Real Points in the Polar Varieties 

Let fx, . . . , fp £ Q [Xx, ■ ■ ■ , X n ] be a reduced regular sequence and let again 
Wo := V(fx, • • • , f p ) be the affine variety defined by fx, ■ ■ ■ , f p ■ Consider 
the real variety So := Wo fl M n and suppose that 

(i) So is nonempty and bounded (and hence compact), 
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(ii) the Jacobian J(/i, . . . , f p )(x) is of maximal rank in all points x of So 
(i.e., So is a smooth subvariety of R n given by the reduced regular 
sequence fi,...,f p ), 

(hi) the variables X±, . . . ,X n are in generic position with respect to the 
polynomials fi,...,f p . 

Further, let C be any connected component of the compact set So , and let 
b := (ai, . . . , a p -i,a p , . . . , a n -i,a n ) G C be a locally maximal point of the 
last coordinate X n in the non-empty compact set C C Sq . Without loss of 
generality we may assume that the upper left (p — 1) -minor m of the Jaco- 
bian J(/i, . . . , / p ) does not vanish in b (by our assumptions there must be a 
(p— 1) -minor of J(/i, . . . , / p ) not vanishing at 6). In any local parametriza- 
tion of So at 6 the variable X n cannot be an independent variable, since 
X n attains a local maximum in b {a n is this local maximum). Hence, with- 
out loss of generality we may assume that the local parametrization of So 
in b has the following form: there exists an open set U C M n_p containing 
the point a := (a p , . . . , a n _i) , and a continuously differentiable function 

ip : U -> R p , ip := (ipi, . . . , ip p -i,ip n ) 

such that 

X\ — tp\ (Xp, . . . , X n _l), • • • , Xp—l — (fip—l (Xp, • • • , X n _l), 
Xn — PniXpj • • • i X n —\) 

holds for any x = (x p , . . . , x„_i) € . With respect to this local parametriza- 
tion, the polynomials fk, 1 < k < p, induce real valued functions of the 
form: 

fk{X p , . . . , X„-i) : = 
fk(fi(X p , . . . ,X n -i), . . . , (f p -i(X p , . . . ,X n -i), 
X p , . . . , X n -i, ip n (X p , . . . , X n -i)). 

For every 1 < k < p, and every p < j < n — 1, one has the identity 
dfk 



dfk_ df^dip^ 
dXj dXj dX x dXj 

in the open set U . 
Considering the (p x p) -matrix 



B :-- 



dfi 
dX 1 



+ 



df P 
dX 1 



dfk 


dcp p -i 


5Xp_x 


dXj 


dfi 


dfi 




dX n 


df P 


df P 



+ 



dX n dXj 



= (1) 



dX, 



p-i 



dX n 
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and observing that B is regular in IA , we obtain from (1) that 



detB(x) 



dpi 



dip 



p-i 



dXj 



(Adj B)(x) 



df 



dX 



df P 



dX 



±(s) 



df P 
dX~ 



(x) 



(2) 



holds for any x €zU (here Adj B denotes the adjoint matrix of the matrix 
B ) . As b is a locally maximal point of X n , we have that 



8X 4 







holds for every p < j < n — 1 . Thus, equation (2) implies 
Bin, 1) (6) |^ (6) + ••• + B(„ >P ) (6) |^ (6) 







(3) 



for every p < j < n — 1 (here for 1 < k < p we denote the entry of the 
adjoint matrix Adj B at the cross point of the k— th column and the last 
row by B(n, k) ). Taking into account the particular form of the matrix B , 
the equation system (3) means that 



det 



df 



dX Y 



dfi 



dX, 



p-i 



(b) 



dh_ 
dX j 



(6) 



df P 



df P 



(b) 



df P 



(6) 



(4) 



axt y ~ J "' ax p _i v ' dXj 

holds for every p < j < n — 1 . Using our notations for the p -minors of the 
Jacobian J(fi, . . . , f p ) , we reinterprete now the equations (4) as 



M p (b) 



M n _i (b) = 0. 



Since m(b) 7^ holds by assumption, Proposition || implies that b belongs 
to the localized determinantal variety (A n _ p ) m . Therefore, we have b S 



WoD(A n _p) m , i.e., the last formal polar variety W n - P contains the point b. 
On the other hand, b is a nonsingular point of Wo and belongs therefore to 
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W n . p = W n . p \ SingW . Thus W n - p is a non-empty set of dimension zero 
that contains the real point b of the arbitrarily chosen connected component 
C of the real variety S . In particular, b G W n _ p n R n C Wj n R n = 5 f 
holds for any 1 < i < n — p . 

These considerations imply the following result: 
Theorem 10 

Let Wo := V(fi, . . . , f p ) be as in Theorem^. If the real variety Sq := 
WoflU" is non-empty, bounded and smooth, and if the variables X%,..., X n 
are in generic position with respect to fi, ■ ■ ■ , f p , then every real formal 
polar variety Si = W{ n M™, 1 < % < n — p, is a non-empty, smooth 
manifold of dimension n — {p + i) and contains at least one representative 
point of each connected component of the real variety Sq . 



3 Real Equation Solving 

The geometric results of Section 2 allow us to design a new efficient procedure 
that finds at least one representative point in each connected component of 
a given smooth, compact, real complete intersection variety. 
This procedure will be formulated in the algorithmic (complexity) model of 
(division-free) arithmetic circuits and networks (arithmetic-boolean circuits) 
over the rational numbers Q . 

Roughly speaking, a division-free arithmetic circuit (3 over Q is an algo- 
rithmic device that supports a step by step evaluation of certain (output) 
polynomials belonging to Q [X\ , . . . , X n ] , say fx,...,f p . Each step of (3 
corresponds either to an input from X\ , . . . , X n , to a constant (circuit pa- 
rameter) from Q or to an arithmetic operation (addition/subtraction or 
multiplication). We represent the circuit (3 by a labelled directed acyclic 
graph (dag). The size of this dag measures the sequential time requirements 
of the evaluation of the output polynomials fx, ■ ■ ■ , f p performed by the 
circuit f3 . 

A (division-free) arithmetic network over Q is nothing else but an arithmetic 
circuit that additionally contains decision gates comparing rational values 
or checking their equality, and selector gates depending on these decision 
gates. 

Arithmetic circuits and networks represent non-uniform algorithms, and the 
complexity of executing a single arithmetic operation is always counted at 
unit cost. Nevertheless, by means of well known standard procedures our 
algorithms will always be transposable to the uniform random bit model 
and they will be practically implementable as well. All this can be done in 
the spirit of the general asymptotic complexity bounds stated in Theorem 



11| below. 

Let us also remark that the depth of an arithmetic circuit (or network) mea- 
sures the parallel time of its evaluation, whereas its size allows an alternative 
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interpretation as " number of processors" . In this context we would like to 
emphasize the particular importance of counting only nonscalar arithmetic 
operations (i.e., only essential multiplications), taking Q -linear operations 
(in particular, additions/subtractions) for cost-free. This leads to the no- 
tion of nonscalar size and depth of a given arithmetic circuit or network f5 . 
It can be easily seen that the nonscalar size determines essentially the total 
size of (3 (which takes into account all operations) and that the nonscalar 
depth dominates the logarithms of degree and height of the intermediate 
results of j3 . 

An arithmetic circuit (or network) becomes a sequential algorithm when we 
play a so-called pebble game on it. By means of pebble games we are able to 
introduce a natural space measure in our algorithmic model and along with 
this, a new, more subtle sequential time measure. If we play a pebble game 
on a given arithmetic circuit, we obtain a so-called straight line program 
(sip). In the same way we obtain a computation tree from a given arithmetic 
network. For more details on our complexity model we refer to p| , |25[ |, 
[ 26 1 , |f45| , |)^ ], [ [3§| 1 and especially to |33|] (where also the implementation 



aspect is treated). 

In the next Theorem |ll| we are going to consider families of polynomials 
fi, . . . , f p belonging to Q [Xi, . . . , X n ] , for which we arrange the following 
assumptions and notations: 

(i) fx, . . . , f p form a regular sequence in Q [X±, . . . , X n ] , 

(ii) for every 1 < k < p the ideal (/i, . . . , fk) generated by fi,...,fk in 
Q [X±, . . . , X n ] is radical and defines a subvariety of (D™ of dimension 
n — k that we denote by := V(/i> • • ) fk) ■ 

(iii) the variables Xi, . . . ,X n are in generic position with respect to the 
polynomials fx,...,f p . 

Let W := {x G C n |/i(x) = ••• = f p (x) = 0} and denote by SingW 
the singular locus of Wq . For 1 < i < n — p let Wi be the i -th formal 
polar variety associated with Wq and the variables X p+ i, . . . ,X n , and let 
Wi :=Wi \ SingWo be the i-th polar variety of Wo in the usual sense (see 
Section 2 for precise definitions). Further, for 1 < k < p we shall write 
% :=V k \ SingWo . We call 

5 := max{max{deg 1 1 < k < p}, maxjdeg Wi\l < i < n — p}} 

the degree (of the real interpretation) of the polynomial equation system 
fi, . . . , f p . Finally, let us make the following assumption: 

(iv) the specialized Jacobian J(/i, . . . , fp)(x) has maximal rank in any 
point x of So := W n M n = {x € M n |/i(x) = • • • = f p (x) = 0} and 
St) is a bounded semialgebraic set (hence, So is empty or a smooth, 
compact real manifold of dimension n — p\ see Section 2 for details). 
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Theorem 11 

Let n,p, d, 5, L and I be natural numbers with d > 2 and p < n . There ex- 
ists an arithmetic network J\f over Q of size ( p "]_) L(nd6)°^' and nonscalar 
depth 0(n(lognd + £) log 5) with the following property: Let fi,---,f p be 
a family of n -variate polynomials of a degree at most d and assume that 
fx , . . . , f p are given by a division-free arithmetic circuit /3 in Q [X± , . . . , X n ] 
of size L and nonscalar depth I . Suppose that the polynomials fi,---,f p 
satisfy the conditions (i), (ii), (in) and (iv) above. Further, suppose that 
the degree of the real interpretation of the polynomial system fi, ■ ■ ■ ,f p is 
bounded by 5 (let us now freely use the notations just introduced before). 

The algorithm represented by the arithmetic network Af starts from the 
circuit (3 as input and decides first whether the complex variety W n - P is 
empty. If this is not the case, then W n - P is a zero-dimensional complex 
variety and the network N produces an arithmetic circuit in Q of asymp- 
totically the same size and nonscalar depth as M , which represents the 
coefficients of n + 1 univariate polynomials q,pi, ■ ■ ■ ,p n € Q [X n ] satisfying 
the following conditions: 

degq = # W n -p, 

max{degpfc|l < k < n} < degq, 

W n -p = {(px(u),... ,p n (u))\u G <C,q{u) = 0}. 

Moreover, the algorithm represented by the arithmetic network J\f decides 
whether the set W n - P fl H™ is empty. In this case we conclude Sq = 
Wo fl R n = . Otherwise, the network M produces at most # W n - P < 5 
sign sequences belonging to the set {—1, 0, 1} such that these sign sequences 
encode the real zeros of the polynomial q "a la Thorn" (fii!\l). In this way, 
namely by means of the Thorn encoding of the real zeros of q and by means 
of the polynomials p\ , . . . , p n , the arithmetic network M describes the finite, 
non-empty set 

W n - P n R n = {(pi(«), . . . ,p n {u))\u e TR,q{u) = 0}, 

which contains at least one representative point for each connected compo- 
nent of the real variety 

S = {xGM n \f 1 (x) = --- = f p (x) = 0}. 

Proof 

We shall freely use the notations of Section 2. Any selection of indices 
1 < i\ < ■ ■ ■ < i p < n and 1 < j, k < p determines a p- minor M(i\, . . . , i p ) 
and a (p — 1) - minor m(i\ , . . . ,i p ;j,k) of the Jacobian J{f\ , . . . , f p ) in the 
following way: M(i\, . . . ,i p ) is the determinant of the (p x p)- submatrix 
of J(/i, . . . , f p ) with columns i\, . . . , i p , and m(ii, . . . , i p ; j, k) is the deter- 
minant of the matrix obtained from the former one deleting the row number 
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j and the column number if, . There are p 2 Q) such possible selections. Let 
us fix one of them, say %\ := 1, . . . , i p := p; j := p,k := p. Then, using the 
notations of Section 2, we have m(ii, . . . , i p ;j, k) = m, M{i\, . .. ,i p ) = M p . 
Let us abbreviate g := mM p . From our assumptions on fi,...,f p and 



Theorem |8] and Theorem 10 of Section 2 we deduce the following facts: For 
any 1 < i < n — p the polynomials fx, ■ ■ ■ ,f p , M p , . . . , M p+ j_i have degree 
at most pd. They generate the trivial ideal or form a regular sequence in 
the localized Q -algebra Q [X%, . . . , X n ] g . In either case the ideal generated 
by fx, . . . , f p , M p , . . . , M p _|_j_i in Q [X±, . . . , X n ] g is radical and defines a 
complex variety that is empty or of degree 



deg(Wi \ V(g)) < deg(Wi \ SingW ) = degWi < 5. 

Moreover, by assumption, the polynomials fi,---,f p form a regular se- 
quence in Q [Xi, . . . ,X n ] g and for each 1 < k < p the ideal generated by 
fl, ■ ■ ■ , fk in Q [X±, . . . , X n ] g is radical and defines a complex variety of 
degree 

deg(V k \V(g)) <degV k <5. 

One sees easily that the polynomials f%, . . . , f p , M p , . . . , M n _i and g can 
be evaluated by a division-free arithmetic circuit of size 0(L + n 5 ) and 
nonscalar depth 0(logn + £) . Applying now, for each 1 < i < n — p, 
the algorithm underlying pC[] , Proposition 18 in its rational version |31], 
Theorem 19 to the system 

/i = 0, . . . , f p = 0, M p = 0, . . . , M p+i _i = 0, g + 

we are able to check whether the particular system 

h = 0, . . . , f p = 0, M p = 0, . . . , M n _x = 0, g ^ 

has a solution in (D 71 . If this is the case, then this system defines a zero- 
dimensional algebraic set, namely W n - p \V(g) , and the algorithm produces 
an arithmetic circuit 7 in Q that represents the coefficients of n+1 univari- 
ate polynomials q,p\, . . . ,p n € Q [X n ] satisfying the following conditions: 

degq = #(W n . p \V(g)), 

maxjdeg pk\l < k < n} < degq, 

W n - P \ V(g) = {(pi(«), . . . ,pn(u))\u G C, q{u) = 0}. 

The algorithm is represented by an arithmetic network of size L(ndS)°^ 
and nonscalar depth 0(n(log nd + t) log 5) , and the circuit 7 has asymptot- 
ically the same size and nonscalar depth. Running this procedure for each 
selection 1 < i\ < ■ ■ ■ < i p < n and 1 < j,k < p we obtain an arithmetic 
network Mq of size p QL(nd5) ^ = ( ")L{nd5)°W and nonscalar depth 
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0(n(log nd + £) log 5) , which decides whether W n - P = W n - P \ SingWo is 
empty. Suppose that this is not the case. Then Mo describes locally the va- 
riety W n - P , which is now zero-dimensional. Each local description of W n - P 
contains an arithmetic circuit representation of the coefficients of the mini- 
mal polynomial of the variable X n with respect to the corresponding local 
piece of W n - P . Moreover, one easily obtains the same type of information 
for any linear form Xi+X n and any variable Xi with 1 < % < n. One multi- 
plies now all minimal polynomials of the variable X n obtained in this way. 
Making this product squarefree (see e.g [45], Lemma 12) one obtains the 
polynomial q of the statement of Theorem 1 1 . Doing the same thing for the 
minimal polynomials of each linear form Xi + X n and each variable Xi with 
1 < i < n, yields by means of 45], Lemma 26, the polynomials pi, . . . ,p n 



of the statement of Theorem 11. All this can be done by means of an arith- 
metic network Mi , which extends Mo and has asympotically the same size 
and nonscalar depth. The desired arithmetic network M is now obtained 
from Mi in the same way as in the proof JTJ, Theorem 8, namely as follows: 
applying the main algorithm of |J] or [61] and adding suitable comparison 
gates for rational numbers, we extend Mi to a new arithmetical network M 
of asymptotically the same size and depth, such that M decides whether 
the univariate polynomial q has a real zero. If this is the case, the network 
M enumerates the existing real zeros of q , encoding them " a la Thorn" 
(|]l8||). If q has no real zero we conclude So = 0. Otherwise, the network 
M encodes all real zeros of q by means of #W n - p < 5 sign sequences be- 
longing to the set {—1, 0, 1} . This encoding and the polynomials pi, . . . ,p n 
describe now the set W n - P n R" = {(pi(u), . . . ,p n (u))\u G M,q(u) = 0} 
that contains a representative point for each connected component of Sq . 
□ 



Remark 12 

(i) Using the refined algorithmic techniques of p^ / or J3^ / it is not too 
difficult to see that for inputs fi, . . . , f p represented by straight-line 
programs of length T and space S the arithmetic network M can be 
converted into an algebraic computation tree which solves the algorith- 
mic problem of Theorem in time 0((Tdn 2 + n 5 )5 3 log 3 5 log 2 log 6) 
and space 0(Sdn5 2 ) . 

(ii) The smooth, compact hypersurface case (with p := 1 ) of Theorem |J 
corresponds exactly to m, Theorem 8. 

(Hi) Let J(fi, . . . , f p ) T denote the transposed matrix of the Jacobian 
J(fi, . . . , f p ) of the polynomials fi, . . . , f p in the statement of Theo- 
rem [II] and let 



V:=detJ(fi,...,f p )J(fi,...,f p ) T . 
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From the well-known Cauchy-Binet formula one deduces easily that, 
with the notations of Section 2, the identity 

V= Yj det 2 M(h,...,i p ) 

l<il<---<i p <n 

holds. Replacing now, in the statement and the proof of Theorem [TJ 
for 1 < i < n — p , the polar variety Wi by Wi := Wj \ V(V) and the 
parameter 5 by 

5 := max{max{deg < k < p}, max{deg Wi\l < i < n — p}} 

one obtains a somewhat improved complexity result, since 5 < 5 holds. 

Let us now suppose that the polynomials fi, ■ ■ ■ , f p € Q [Xi, . . . , X n ] sat- 
isfy the conditions (i), (ii), (iii), (iv) above. Unfortunately, the complexity 
parameter 5 of Theore m jlTl is strongly related to the complex degrees of the 
polar varieties W u ..., W„- p of W = {x 6 C n |/i(x) = • • • = f p (x) = 0} 
and not to their real degrees. Under some additional algorithmic assump- 
tions, which we are going to explain below, we may replace the complexity 
parameter 5 by a smaller one that measures only the real degrees of the 
polar varieties W\, ■ ■ ■ , W n - P . We shall call this new complexity parameter 
the real degree of the equation system f\ , . . . , f p and denote it by 5* . 
Let 1 < k < p and let us consider the decomposition of the intermedi- 
ate variety V/% into irreducible components with respect to the Q -Zariski 
topology of <C n say Vfc = C\ U • • • U C s . We call an irreducible component 
C r , 1 < r < s , real if C r n H n contains a smooth point of C r . The union 
of all real irreducible components of is called the real part of and 
denoted by U fc * . We call deg U fc * the real degree of the intermediate variety 
Vk . Similarly, we introduce for every 1 < i < n — p the real part W* of 
the polar variety Wi and its real degree deg W* . Finally, we define the real 
degree of the equation system fi,...,f p as 

(5* := max{max{deg V£\l < k < p},max{degPU*|l <i<n — p}}. 



Now, we are going to restate the main outcome of Theorem 11 in terms of 
the new complexity parameter 5* . For this purpose we have to include the 
following two subroutines in our algorithmic model: 

• the first subroutine we need is a factorization algorithm for univari- 
ate polynomials over Q . In the bit complexity model the problem 
of factorizing univariate polynomials over Q is known to be of poly- 
nomial time complexity |pl|| , whereas in the arithmetic model we are 
considering here this question is more intricate p7p . In the extended 
complexity model we are going to consider here, the cost of factorizing 
a univariate polynomial of degree D over Q (given by its coefficients) 
is accounted as . 
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• the second subroutine allows us to discard non-real irreducible com- 
ponents of the occurring complex polar varieties. This second sub- 
routine starts from a straight-line program for a single polynomial in 
. . . , X n ] as input and decides whether this polynomial has a real 
zero (however, without actually exhibiting it if there is one). Again 
this subroutine is taken into account at polynomial cost. 



We call an arithmetic network over Q extended if it contains extra nodes 
corresponding to the first and second subroutine. 

Modifying our algorithmic model in this way, we are able to formulate the 
following complexity result, which generalizes Q, Theorem 12 and improves 



the complexity outcome of our previous Theorem 1 1 



Remark 13 

Let n,p,d,5* , L and t be natural numbers with d > 2 and p < n. There 
exists an extended arithmetic network AT* over Q of size ( " 1 )L(nc?(5*) ( 1 - ) 
with the following property: Let fi, ■ ■ ■ , f p be a family of n-variate poly- 
nomials of a degree at most d and assume that fx , . . . , f p are given by a 
division-free arithmetic circuit (3 in Q \X\ , . . . , X n ] of size L . Suppose 
that the polynomials fi, . . . , f p satisfy the conditions (i), (ii), (in), and (iv) 



contained in the formulation of Theorem Let us now freely use the 
notations introduced in the present section. Assume that the real variety 
So = {x £ H n |/i(x) = • • • = f p (x) = 0} is not empty and that the real 
degree of the polynomial system fi,...,f p is bounded by 5* . The algo- 
rithm represented by the arithmetic network N* starts from the circuit (3 
as input and decides first whether the complex variety W*_ p is empty. If 
this is not the case, then W*_ p is a zero-dimensional complex variety and 
the network Af* produces an arithmetic circuit in Q of asymptotically 
the same size as Af* , which represents the coefficients of n + 1 univariate 
polynomials q* ,p\, . . . ,p* n € Q [X„] satisfying the conditions 

degq* = #W:_ p , 

max{degp|,|l < k < n} < degg*, 

W:_ p = M(u), . . . ,p* n (u))\u G C, q*(u) = 0}. 

Each over Q irreducible component of the complex variety W*_ p contains at 
least one real point characterized by an irreducible factor of the polynomial 
q* . The algorithm represented by the network Af* returns all these points 
in a codification "a la Thorn". Moreover, the non-empty set W*_ p D H n 
contains at least one representative point for each connected component of 
the real variety Sq . 



The proof of this remark is a straight -forward adaptation of the arguments 
of the proof of Theorem 12 (which treats only the hypersurface case 
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with p := 1) to the arguments of Theorem [0] above. Therefore, we omit 
this proof. 

Let us finally observe that the practical relevance of the complexity out- 
come of Remark 6 is highly hypothetical, because it depends on the strong 
assumption that extended arithmetical networks are realizable by perfor- 
mant, programmable algorithms. Nevertheless, by means of Remark 6, we 
wish to underline the importance of the search for efficient procedures that 
realize the first and second subroutine introduced as extra nodes in our 
complexity model of extended arithmetic networks. 
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